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Radiative Penguin Decays from BaBAR 
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We summarize the latest BABAR results on B ^ and B p(cS)y- 



1 Introduction 

Electroweak penguin decays provide a promising hunting 
ground for Physics beyond the Standard Model (SM). The 
decay B Xsj, which proceeds through an electromag- 
netic penguin loop, already provides stringent constraints 
on the supersymmetric (SUSY) parameter space [ 1 1. The 
present data samples of ~ 1 x IQ^BB events allow to ex- 
plore radiative penguin decays with branching fractions of 
the order of 10"^ or less. In this brief report we discuss 
a study of B — > K'^*^i^{^ decay modes and a search for 
B p(to)y decays. 

2 Study of 5 ^ Krr and B K*rr 

The decays B KC^r and B K*e{- proceed through 
an electromagnetic penguin loop, a Z° penguin loop or a 
weak box diagram as shown in Figure[n In the framework 
of the operator product expansion (OPE) the decay rate 
is factorized into perturbatively calculable short-distance 
contributions that are parameterized by scale-dependent 
Wilson coefficients and non-perturbative long-distance ef- 
fects that are represented by local four-quark operators. 
Operator mixing occurring in next-to-leading order pertur- 
bation theory leads to three effective scale-dependent (/i^) 
Wilson coefficients, C'^^^ (jib),Cg-^^(fj.h), and C^^-j (jib) that 
are each sensitive to New Physics contributions. Exam- 
ples for non-SM penguin loops are depicted in Figure |2l 
In the Standard Model the branching fractions are pre- 
dicted to be within the following ranges, S(B — > Kt^C^) - 
(0.23 - 0.97) xlO-'',S(B ^ K*n^^r) = (0.81-2.64)xl0-'' 
and S{B K*e+e-) = (1.09 - 3.0) x 10"^ [ 2|. In su- 
persymmetric models, for example, the branching fractions 
my be enhanced by more than a factor of two [ 1 1 . 

BABMi has analyzed eight final states where a K^, K^^,K*^^ 
or K*"^ recoils against a ^i^iJT or e^e^ pair, us- 
ing an integrated luminosity of 77.8 fb"' that corre- 
sponds to (84.4 + 0.9) x lO*" BB events. The dis- 
criminating variables are beam energy-substituted mass 

niEs = ■\l(Ei,emn'^^ ~ (Pb>^ ^"^^ *^ energy difference 
AE* - Et - El , where pn, E*„ and El denote the 

D beam ^ d beam 

B-momentum, B-energy and beam energy in the center- 
of-mass (CM) frame, respectively. The AE* - niEs plane 
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Figure 1. Lowest-order diagrams for B — > A^***^+f" in SM. 
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Figure 2. Penguin loop diagrams involving a W boson, a charged 
Higgs boson, a chargino and a neutralino, respectively. 

is divided into three regions, a signal region (±3cr boxes 
around the signal), the fit region (niEs > 5.2 GeV, 
I AE* |< 250 Me V), and a large side band (m^s > 5.0 GeV, 
I AE* |< 500 MeV). Specific selection criteria are used to 
suppress individual backgrounds. A Fisher discriminant [ 
^5 1 is used to eliminate the back-to-back continuum back- 
ground. It is based on event shape variables such as the 
thrust angle between the daughter particles of the B meson 
candidate and that of the remaining particles in the event 
(cos 6*7-), the B decay angle in the Y(45) rest frame be- 
tween the B candidate and the beam axis (cos 6*^), and the 
ratio of second-to-zero th Fox- Wolfram moments (R2) [ 3|, 
as well as the invariant mass of the K - ( system mu- To 
discriminate against combinatorial BB background a likeli- 
hood function is used that combines the missing energy in 
the event {E,„jss), with the dilepton vertex probability, the 
significance of the dilepton separation along the beam di- 
rection and cos ff'g. To reject events from B —> J/tfrK^*^ and 
B t//(2S)K^*> decays with J/ip(i//(2S) t^T that have 
the same event topologies in a restricted mtc mass region 
as signal events, the shaded regions in the AE* - mtc plane 
shown in Figure|3are vetoed. The inclined bands provide 
an efficient rejection of J/ifrK'^*^ events in the fit region, in 
which one or both leptons radiated a photon. 

The selection criteria are optimized on simulated data as 
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Figure 3. Charmonium veto (hatched regions) in the A£* - m^+f- 
plane for B — > A'**'e+e" and B — » K'-'''\i*\r . The dots repre- 
sent simulations for B Jl^{-^ er)K^''> and B ^ iJ/{2S){^ 

rr)A:<*>. 



well as data sidebands. The simulated data are used to de- 
termine efficiencies and estimate the contribution of back- 
grounds that may peak in the AZs'-rngs signal box. Specif- 
ically, the decay channels B — > K*y, B J/iffK^*\ 
B -> >^(2S)K^*\ B K^*'>7r°, B K^*^t], B Knn, 
B KKn, B KKK, and B ^ Dk have been stud- 
ied. The Monte Carlo results are checked with data control 
samples, including exclusive and inclusive charmonium 
decays, B — > Dii decays, data sidebands and K^*''e'^if sam- 
ples. For the combined exclusive charmonium modes we 
achieve an excellent agreement between data and Monte 
Carlo for the ratio of event yields of 1 .013 ± 0.018. Further 
details of the event selection are discussed in [ 6|. 

In each of the four K{^(^ final states, a signal is extracted 
from a two-dimensional fit to the iues - AZi* plane. Sim- 
ilarly, a signal yield for the four K*(*(^ final states is ob- 
tained from a three-dimensional fit to the niEs - Ais* - mj^^j 
space. The signal shapes are obtained from Monte Carlo 
samples with fine-tuning on the exclusive charmonium 
modes. To account for the asymmetric AZs* and Mes dis- 
tributions that result from radiation effects of the final-state 
leptons and correlations in the two variables a product of 
Crystal Ball functions is used [ Tj. The combinatorial 
backgrounds are parameterized with ARGUS functions [ 
IHl, where both the normalization and the shape parame- 
ters are left free. Only in B — > K'^e^e^ , a significant yield 
of 16.2;^^ J events is observed. The selection efficiency is 
(19.1 + 1.1)%, yielding a branching fraction of 



S(B+ ^ K^e^e-) = {Q.96^^f^ + 0.06) x 10"^ 



(1) 



The second largest yield of (V.S^^ j) events is found for 
B^ — > K*^e^e^ for which a selection efficiency of (10.6 ± 
0.9)% is achieved. The individual results for each of the 
eight final states are summarized in Tabled 

The Mes and tS.E* projections of the two-dimensional fit 
for the combined B — > KC^£^ final states are shown in 
Figure |3 The coiTesponding distributions of the three- 
dimensional fit for the combined B K*t^i^ channels 
as well as the niKn mass projection are shown in Fig- 
ure |5] In order to combine K'fi^jj'' and K*e^e^ results 
their ratio of branching fractions is assumed to be S{B — > 
K*e^e-)/!B(B Ky^n') = 1.33 [ 1|. The multiplica- 
tive systematic errors that do not affect the significance 
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Figure 4. Projections of the two-dimensional fit of the combined 
KC^C^ data sample for a) iiies after the constraint -0.11 < AE* < 
0.05 GeV and b) A£* after requiring \mEs - triBl < 6.6 MeV/c~. 
The solid and dashed curves show the entire fit and the total back- 
ground contribution, respectively. 
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Figure 5. Projections of the tree-dimensional fit of the combined 
K'H^C" data sample for a) niEs after the constraint -0. 11 < AE* < 
0.05 GeV and 0.817 < rriKn < 0-967 MeV, b) ^E' after requiring 
\mEs - mB\ < 6.6 MeV/c^ and 0.817 < niKn < 0.967 MeV, and 
c) niKn satisfying -0.11 < AE* < 0.05 GeV and \mEs - nis] < 
6.6 MeV/c^. The solid and dashed curves show the entire fit and 
the total background contribution, respectively. 



range from 6% - 12% for Kt^T modes and 7% - 15% for 
K*{^{^ modes. The largest contributions result from the 
K'^ efficiency (10%), the model dependence (4% - 7%), fx 
identification (3.2%), BB likeHhood ratio (1.3% - 5.6%), 
the Fisher discriminant (0.9% - 3.1%), and the tracking ef- 
ficiency for hadrons (1.3% - 3.9%). Contributions from 
K/n identification, e identification, tracking efficiency for 
leptons, Monte Carlo statistics and BB counting are small 
(1% - 2%). The additive systematic errors that affect the 
significance result from the signal yields in the fit, includ- 
ing uncertainties in signal shapes, background shapes and 
the amount of peaking backgrounds. Including systematic 
errors the significance of the combined Kl^i^ sample is 
7. Oct, while that of the combined K*{^{^ sample is 3.0cr. 



For the combined B 



Kl^{ modes we measure a 
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Table 1. Measured event yields, efficiencies and brandl- 
ing fractions for B — > K£*£^ and B — » K*(*C^ modes. 



Mode 


yield 
[events] 


e[%] 


Sx 10^ 




16.2!:-» 
-0 5+'-'' 

^■-'-1.0 
j.u 2 

9 7+4.1 
z,./_2g 
4 4+3.5 


19.1 + 1.1 
8.7 + 0.6 
20.6 + 2.5 
9.2 + 0.7 
13.5 + 0.9 
6.5 + 0.6 
10.2+ 1.4 
5.0 + 0.8 


n q-r +0.28+0.06 
'-'•^"-0.24-0.06 
n 14+0.25+0.02 
^•^^-0.16-0.02 
n no+0.25+0.05 

1 70+0.99+0.22 
^•'"-0.73-0.22 
n no+0.68+0.18 
^•^"-0.54-0.18 
1 10 +0.99+0.26 

^•-'^-1.38-0.33 
A TT+3.4I+0.68 
^-'•'-7,.40-O.fiS 



branching fraction of 



S(B ^ Ke^e-) = (0.68!^ !^ + 0.04) x 10"^ 



(2) 



For the combined B K*{*( modes we determine a 
branching fraction of 



^ K'e^e-) = (1.4!jjj^ + 0.21 x 10"^ 



(3) 



The BMi/R results are consistent with those obtained by 
BELLE and most SM predictions, but S{B Kt^C) 
is almost 2cr higher than a recent next-to-next-to-leading 
order prediction of ^ Kl^r) = (0.35 +0.12) x 10"^ 

[ID. 



3 Search for B p(.co)y 

The decays B —> p{(jS)y are flavor-changing neutral cur- 
rent b d transitions that are mediated by an electro- 
magnetic penguin loop. The decay rates depend on the ef- 
fective Wilson coefficient C'i/^fib) that may be enhanced 
by New Physics contributions as in B — > K*y. They 
are, however, suppressed with respect to B ^ K*y by 
WidlVts?'- Thus, measuring the ratio of branching frac- 
tions S{B — > py)/S(B — > K*-/), allows us to extract 
|Vfrf/V/.5|. In SM the branching fraction for the charged 
mode in next-to-leading order is predicted to lie between 
S(B+ ^ p+y) = (0.85 ± 0.4) x lO"*^ [[El and S(B ^ 
p+y) = (1.58+jJ-^^) X 10"*' [|13|. For the neutral modes the 
branching fractions are a factor of two smaller than that 
for B^ — > p^j due to isospin. By determining the ratio of 
decay rates r(B py)IY{B K*y) uncertainties due to 
scale dependence and a precise knowledge of form factors 
cancel except for SU(3) breaking effects. The largest re- 
maining theoretical uncertainties result from contributions 
of W-annihilation/W-exchange diagrams, which amount to 
~ 20% in 5+ p+y and to ~ 13% in B° ^ p°(w)7 r lHl . 



9.6 fb ' in the continuum 40 MeV below the T(45) peak. 
Challenges in the analysis stem from a huge qq continuum 
background including initial-state radiation, background 
from B K*y and the fact that the p resonance is much 
broader than the K* resonance. The qqy continuum back- 
ground with a hard y from initial-state radiation may have 
a similar event shape as that of the signal which is less 
spherical than a typical BB event. Photon candidates with 
energies of 1.5 GeV < Ey < 3.5 GeV that are inconsis- 
tent with originating from a tt" or 77 decay are combined 
with a p*,p^, or oj candidate, where the vector mesons 
are reconstructed from 7!-^;7r'\ tt^tt^ and three-pion combi- 
nations, respectively. Rejecting charged tracks in the sig- 
nal that are consistent with a kaon, the K/n misidentifica- 
tion is less than 1%. The ;7r;7r (3;7r) invariant mass has to 
lie within a mass window of 520 - 1020 MeV/c^ (759.6 - 
805.6 MeV/c^) and its momentum in the CM frame must 
satisfy 2.3 < pj^ < 2.85 GeV/c (2.4 < p^^ < 2.8 GeV/c). 
A tt" candidate must have a yy invariant mass of 115 < 
niyy < 150 MeV/c^. To improve the momentum resolu- 
tion we perform a kinematic fit with niyj. constrained to the 
nominal tt" mass. 

To reduce the qq continuum background a neural net- 
work is used that is based on event-shape variables 
(cos Othrust, cos 0*^, cos Oheiicity), the Dalitz decay angle for 
cj, the energy flow in 18 cones around photon direction, 
the vertex separation Az, the ratio of second-to-zeroth Fox 
Wolfram moment, R'2, calculated in a frame recoiling the 
photon and the net flavor in the event F 1101 . The neural 
network is trained with Monte Carlo signal events and con- 
tinuum data. The neural network output is cross-checked 
with data using a Z?" D n^ sample that has a similar 
topology as the signal and the off-resonance data sample. 
Further details of the event selection are discussed in [ 1 IJ. 

The AE* - niEs distributions for the final data samples are 
shown in Figure |6l The signal yields are extracted from 
a maximum likelihood fit in the three-dimensional space 
AE* - niEs - mpinitj). The procedure is crosschecked with 
our B -> K*y data samples. For B — > K*'^ {K*^)y the maxi- 
mum likelihood fit yields 343.2 + 21.0(93.1 ± 12.6) events, 
which is in good agreement with the expected yields of 
332 ± 36 (105 + 18) events, respectively. 

The extracted signal yields of 4.8^H events for B p^y. 



6.2^^ T events for B 

— 6.Z 



BAB/U? has searched for B p{(^)y modes using an in- 
tegrated luminosity of 77.8 fb ' on the T(45') peak and 



p^y, and 0.1;^2o events for B — > 
a>y are consistent with background. The efficiencies are 
12.3%, 9.2% and 4.6%, respectively. Including systematic 
errors, which respectively increase from 11.8% to 13.4% 
and 17.3% for the three decay channels, we obtain branch- 
ing fraction upper limits @ 90% CL of S(B" ^ pV) < 
1.2x 10"^ 

S(B+ p^y) < 2.1 X 10-^ andS(B" ^ wy) < 1.0 x 10"^ 
These limits are significantly lower than those of previous 
searches rlT6lrrT71. 

Assuming isospin symmetry the p^y and p"y samples are 
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Figure 6. AE' - ities scatter plots of the fit region for a) — » 
p^y, b) S" — > p^y and c) B" — > a>y candidates. The boxes indicate 
the expected signal regions 



Table 2. Extrapolations of the significance, experimental error on 
the branching fraction and experimental error on |V,f;/V,.tl for the 
combined B — » p(cj)y modes expected for different luminosities. 



Luminosity 


significance 


(o-s/S)„p 


(r(V,JV„) 


100 fb-i 


1.9 


-2.8 


cr 


0.38-0.53 


0.19-0.27 


200 fb-i 


2.7 


-3.9 


cr 


0.28-0.38 


0.14-0.19 


300 fb-i 


3.3 


-4.8 


cr 


0.23-0.31 


0.12-0.15 


400 fb-i 


3.9 


-5.5 


cr 


0.20-0.27 


0.1-0.14 


500 fb-i 


4.3 


-6.2 


cr 


0.18-0.25 


0.09-0.13 


1000 fb-i 


6.0 


-8.7 


0" 


0.14-0.18 


0.07-0.09 



on \V,d/V,s\) result from the present event selection, while 
the large (small) values are obtained by assuming that for 
the same selection efficiency as in the present analysis the 
background is halved. For a luminosity of ~ 500 fb^^ a 
significant measurement of these modes is expected. The 
limiting factor for extracting IV/rf/y,.,! is the theoretical un- 
certainty from S U(3) breaking and the size of the annihi- 
lation diagram. The latter uncertainty can be removed by 
using only Z?" — > p'^y events, but for the same significance 
a factor of four increase in luminosity is required. 



4 Outlook 



combined yielding !B(B py) < 1.9 x 10"^. Using the 
recent BABAR B — > K*y branching fraction measurement 
[ 1141 this translates into an upper limit on the ratio of 
branching fractions of S(B ^ py)/S(B K*y) < 0.047 
@ 90% CL. To constrain the ratio of |Vfd/V,.5| we use the 
parameterization F fT^ . 



g(g ^ py) 
!B(B K*y 



■ ml/Ml ^ 



I -ml, I Ml) 



(^{\ +AK\. 



(4) 



The parameter ^ represents S U(3) breaking while AR ac- 
counts for the annihilation diagram in B^ —> p^y. Using 
f = 0.76 + 0.1 and AR ^ Q.Q + 0.2 ([[T3, [[151) we ob- 
tain an upper limit of \V,dlV,s\ < 0.34 @ 90% CL. This 
is still larger than the limit of \V,dlV,s\ < 0.23 @ 95% CL 
which is derived from BsBs and BdBd mixing results for 
Am/j, = 0.503 + 0.006 ps"', Am^, > 14.4 ps"' @ 95% CL 
and^= 1.24[[TH|. 

The present upper limits are approaching the theoretical 
predictions. Assuming a branching fraction of 
p^y) = 1 X 10"^ and ^ pV = S{B or/ ^ 

jS(B^ — > p^y) we estimate a signal significance and 
experimental errors for different luminosities as listed 
in Table |2 The small values for the significance (i.e. 
large experimental errors on the branching fraction and 



By 2007 BABAR expects to record an integrated luminosity 
of ~ 500 fb"'. This sample will be sufficient to measure 
the B — > K'-'^i^H^ and B — > p(w)y branching fractions with 
reasonable precision. Due to the theoretical uncertainties, 
however, tests of the Standard Model will be rather limited. 
For example, at a luminosity of 1 ab"' the presently quoted 
theoretical uncertainties for |y,rf/y,j| will be larger than the 
extrapolated experimental errors. 

Observables that are barely affected by theoretical uncer- 
tainties and, therefore, provide excellent tests of the Stan- 
dard Model are the lepton forward-backward asymmetry 
as a function of mc( measured in B — > K^*^{^(' modes and 
direct CP violation in B ^ py channels. In SM the lepton 
forward-backward asymmetry in the B rest frame for dilep- 
ton masses below the J/i// has a characteristic shape, cross- 
ing zero at a specific dilepton mass [[Q. The zero point 
is predicted in SM with small uncertainties. With limited 
statistics we will just determine the zero point, while with 
high statistics we will be able to measure the entire distri- 
bution. Deviations from the SM shape will hint to New 
Physics. In SM CP asymmetries in B ^ py could be as 
large as 12%, but may be modified considerably in models 
with minimal flavor violation [121. Precise measurements 
of the shape of the lepton forward-backward asymmetry 
and direct CP violation in B — > py, however, require data 
samples that are several tens of ab"' . 
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